Instabilities in neutrino systems induced by interactions with scalars 
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If there are scalar particles of small or moderate mass coupled very weakly to Dirac neutrinos, 
in a minimal way, then neutrino-anti-neutrino clouds of sufficient number density can experience 
an instability in which helicities are suddenly reversed. The predicted collective evolution is many 
orders of magnitude faster than given by cross-section calculations. The instabilities are the analogue 
of the "flavor-angle" instabilities (enabled by the Z exchange force) that may drive very rapid flavor 
exchange among the neutrinos that emerge from a supernova. These exchanges do require a tiny 
seed in addition to the scalar couplings, but the transition time is proportional to the negative of the 
logarithm of the seed strength, so that the size of this parameter is comparatively unimportant. For 
our actual estimates we use a tiny non-conservation of leptons; an alternative would be a neutrino 
magnetic moment in a small magnetic field. The possibility of a quantum fluctuation as a seed 
is also discussed. Operating in the mode of putting limits on the coupling constant of the scalar 
field, for the most minimal coupling scheme, with independent couplings to all three v, we find a 
rough limit on the dimensionless coupling constant for a neutrino-flavor independent coupling of 
G < 10~ 10 , to avoid the effective number of light neutrinos in the early universe being essentially 
six. If, on the other hand, we wish to fine-tune the model to give a more modest excess (over three) 
in the effective neutrino number, as may be needed according to recent WMAP analyses, it is easy 
to do so. 

PACS numbers: 13.15.+g 



1. INTRODUCTION 

In the last few years there has been an accumulating 
literature focused on the possibility that the neutral cur- 
rent v — v interaction can induce rapid collective flavor 
exchanges of physical interest in situations in which there 
are initial flavor-energy correlations or flavor-angle corre- 
lations (H-[I3l- Up to now the application has been to the 
supernova neutrino pulse, where, apart from the dream 
that a galactic event will bring us information on neu- 
trino properties directly from terrestrial detection, the 
flavor-energy correlations may matter to the explosion 
dynamics or to the possible R-process nucleosynthesis on 
the periphery of the supernova. 

The purpose of the present paper is to show that the 
same kind of effect can enter in models in which neutri- 
nos interact with scalar particles. Such interactions have 
been the subject of much speculation [l5j]-[25j]. We be- 
lieve that the exploration of different models of such an 
interaction will lead to a whole zoo of exotic phenomena, 
but we focus on the simplest example in this preliminary 
study, and we shall consider only one flavor of Dirac neu- 
trino coupled to a single scalar field. The effects that we 
find are of most interest with respect to the early uni- 
verse, in the period immediately prior to light element 
nucleosynthesis. In this system we have a thermal distri- 
bution of left-handed neutrinos, and their anti-particles, 
in nearly equal number. In standard theory a negligible 
number of right handed light neutrinos have been pro- 
duced up to that early time, through a combination of 
collisions and the tiny neutrino mass term. 

The collective effect that can come into play through 
the medium of the scalar particle exchange is the sud- 



den and simultaneous creation of right-handed i/'s and 
left handed P's, keeping the lepton number equal to zero, 
and with a rate orders of magnitude greater than one 
would predict from a perturbative cross-section. The 
constraint on the model would thus be the need to avoid 
excessive production of right-handed v's and their an- 
tiparticles prior to neutrino decoupling in the big bang; 
that is, to avoid increase of the effective number of neu- 
trino species to above the limits set by cosmology. In 
this system the collective effects of the v — v force arising 
from Zq exchange vanish, owing to the cancellation be- 
tween the particle and anti-particle terms in the effective 
potential. 

As in the supernova case, the surrounding medium in 
our early-universe environment is filled with other par- 
ticles than neutrinos. We shall consider only cases in 
which collision times for the neutrinos to scatter on these 
other particles are long compared to the instability times. 
Then the only effect of the other particles is through 
the constant background neutrino potential coming from 
from the neutral current interactions with electrons and 
positrons (the interactions with nucleons being flavor in- 
dependent and therefore irrelevant). In the early universe 
the number of electrons so nearly equals the number of 
positrons that this effect is negligible. 

We first take a rather long detour through the exist- 
ing lore of collective instabilities that derive from the Zq 
force. It is not our intention to write a review article, but 
we do need to discuss in detail what is the same and what 
is different in the scalar-mediated case. Furthermore, we 
wish to underscore two areas in which one might ques- 
tion the assumptions of the whole enterprise, including 
the work in refs. (H-[3]: 1) the replacement of the in- 
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teraction Hamiltonian by a "forward" Hamiltonian; 2) 
after this replacement, and the extraction of equations of 
motion, a "mean field" assumption. In the case of poten- 
tial issue #1, the issue is, in a pragmatic way, de-fused 
(but not put to rest) by staying in domains in which the 
interesting physics takes place over distances much less 
than a mean-free path as calculated from cross-sections. 
In the case of the mean-field assumption, we address the 
question directly from simulations with a finite number 
of neutrinos in a box. Assemblages of a few hundred neu- 
trinos are computationally accessible on a desktop, and 
with this number we can support the mean-field method 
in some of the simpler coupling schemes. We return to 
these issues in section 4. 



2. FORWARD HAMILTONIAN FROM Z 
EXCHANGE. 

In all cases we begin with an effective interaction 
Hamiltonian, Hi, arising from an exchange of a boson. 
Then if we are operating in , say, the MeV region of v en- 
ergies the v — v interaction mediated by Zq exchange can 
be taken as a local four-Fermi interaction. In the case 
of our scalar exchanges we do want to cover the possibil- 
ity of exchanged masses in the region below lMeV and 
thus the propagator for the scalar particle contributes 
momentum dependence. 

The key to all of the literature cited in refs.[l|-[3] is 
the replacement of Hj by a "forward" Hamiltonian. In 
this reduction we expand each field in plane wave modes; 
each term now is characterized by a momentum sequence, 
Pi) P2 —> P3- P4- The forward interaction Hamiltonian is 
defined by omitting from the sum all terms except those 
in which pi = p 3 ; p 2 = P4 or pi = p 4 ; p 2 = P3 Then in 
discussing the evolution of some initial state under the 
influence of Hj we can ignore the kinetic term altogether. 
The only physical effects that ensue are now the exchange 
of flavors among momentum states (or we could equally 
say: the exchange of momenta among flavor states) in 
the case of the Zq model, or the alteration of helicities, 
with momenta unchanged, for the case of the scalar inter- 
actions. In dynamical calculations the free Hamiltonian 
Hq, is now irrelevant; we deal entirely in a subspace of 
states of equal unperturbed energy. 

These "forward" interaction Hamiltonians can now 
be completely expressed in terms of bilinear forms in 
neutrino creation or annihilation operators, a|(p) %-(p), 
where {i,j} indices are of flavor for the Zq case, and of 
helicity for the scalar case, 



Pij(p) = ai(p) f aj(p) 
Pi,j(P) = aj(p) f a 4 (p), 

with the commutation relations, 

[ftj(p),ft,i(p')] = [kiPkj(p) - <Wm(p)]<W 



(i) 



[Pij (p), Pk,i (p')j = [-$i,iPk,j{p) + <5j,fcft,i(p)]<W' ■ ( 2 ) 
In terms of these operators the forward Hamiltonian is, 



H vv {p) 



V2G f 



2>J 1 



COS 



[Pi,j(p) - Pi,j(P)) - Pj.M) 

+{pi,i{v) - PiAp)) (pjj(q) - Pij(q) 



(3) 



where the indices i,j range over the values e,\x,T. 

We divide the momentum space into N B bins, each of 
them containing N v /Nb neutrinos, and distinguished by 
the superscript (a). With few exceptions the instabil- 
ity literature is based on two-flavor models i = {1,2}. 
In this case it is convenient to use angular momentum- 
like collective variables , , , formed from the 
p{p)i,j by summing over the region of momentum space, 
such that p C a, at the same time dividing out a factor 
N V /N B . 
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n b n: 



pCa 

- 1 ^[5t(p)6(p) + a t (p)a(p)] (4) 

pCa 



where a(p), 6(p) are the annihilation operators for the 
respective flavor states. We use a similar construction for 
anti-particle operators; but at this point we specialize to 
the case where none are present, for simplicity in sketch- 
ing out methods. The anti-particles are easily reinstated 
for actual applications. The commutation rules for these 
collective operators are 



[S£\SW]=N B N- I 8 a , p si a) 
[S { + a \sf} = -2N B N- 1 S a jSi a) 



(5) 



etc.; within each beam, as for the Pauli matrices a± — 
(cxi ± z<72)/2 and 03, except for the factor of NbN^ 1 . 
The Hamiltonian in terms of these variables for the case 
of the neutral current interactions is then, 



H = V2G F n u N v N B x ^ [a a ,pS ( * ] S m + b a ^S { 3 a) S : 



(0) 
3 



a,0 



+c a , S^S^ +H.C. 



(6) 



where n v is the neutrino number density and the matrices 
a and b are determined by the geometry of subdivisions 
into beams, taking into account the angular factors (1 — 



COS 6-r 



In the equations of motion that follow from 



© and ©, generically idS/dt = [H, S], the factor of N v 
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from the Hamiltonian cancels the factor of N v 1 from the 
commutation rules and we obtain, e.g., 



an algebra 



dt 



-b a ,,(Sfs^+S^Sf] 



(7) 



The literature cited in refs.(l|-(ll| is based on these 
equations, but with a further sweeping assumption, 
namely that on the right hand side we can replace the 
expectation values of products of 5* operators with the 
products of the individual expectation values, e. g., 



(8) 



generally characterized as the mean- field limit. 

The diagonal elements of the matrices a and b are zero, 
due to the (1 — cos (9) factor in the original Hamiltonian. 
Therefore there are no operator-order ambiguities on the 
right hand side of (J7J. 

Clearly, in the mean field limit when we take pure fla- 
vor states as an initial condition we have (S"(t = 0)) = 0. 
Then the equations ([7]) and their coupled counterparts 
for the operators S% say that nothing whatever happens. 
But if we took an initial values (5"(0)) 7^ 0, or added to 
H an ordinary neutrino oscillation term, then, according 
to the literature a bewildering number of things can hap- 
pen that depend on the interplay of the two parts of the 
Hamiltonian. Beginning in sec. 4 we have studied some 
aspects of these effects perhaps more systematically than 
does the previous literature. 



3. SCALAR PARTICLE EXCHANGE MODEL 

We take the coupling to a scalar field <j> to the Dirac 
neutrinos to be Hi — gipipcj). Now the effective Hamil- 
tonian for the neutrino-neutrino interaction consists of 
a contribution from a scalar particle exchange and one 
from virtual annihilation. 

We shall consider the scalar mass to be of the order of 
the neutrino energy scale or less. The p, q dependence in 
this case is more complex than in the Z exchange case. 
For a each momentum p there are now four states, neu- 
trino and anti-neutrino each with either normal helicity 
(- for i>, + for v) or opposite helicity. For the massless 
neutrino case, the coupling to the scalar connects the 
"normal" state only to the "opposite" helicity. However 
in constructing the forward Hamiltonian we shall need 
all 16 of the bilinear products of a creation and an an- 
nihilation operator, which we designate p a . We use a 
notation in which a(p) annihilates the v state with (nor- 
mal) negative helicity and a(p) the v state of (normal) 
positive helicity, with &(p),6(p) annihilating the states 
of opposite helicity. The bilinears, a^a, aJb, a% etc. have 
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Mp), Pp{v')] = h l,pP-y(P) S p,p' ( 9 ) 
7=1 



The forward effective Hamiltonian for exchange of the 
scalar particle of mass m, written in terms of these oper- 
ators, is, 

ffca ^EE©(^(PW(q)[^(M 



(10) 



where 



£i[p,q] = (p-q) 2 +m 2 -(|p|-|q|) 2 , (11) 



and 



£> 2 [p,q] = (p + q) 2 



|q|) 5 



(12) 



Here the coefficients £ and 77 are independent of the 
directions of p, q, as long as we use helicity states to 
define the p Q (p). The D\ term comes from the exchange 
of the scalar between any pair of v — v pair, D — v pair or 
v— P pair. 1 The D2 term comes from virtual annihilation 
of a v — v pair. 

In contrast to the case of the Zq exchange terms of 
Q that drove instabilities only in the case of nonistropic 
momentum distributions, the form (|10j) will drive insta- 
bilities in perfectly isotropic systems of neutrinos and 
anti-neutrinos. These instabilities are the main point of 
the present paper. 

We can then introduce ab initio angular averages of the 
operators in (TTU|) . However, also in contrast to the model 
of ([3]) , we wish to consider a range of scalar particle mass 
m that ranges from very small, in comparison to the neu- 
trino momenta, to equal or larger. Thus, while lacking 
the cos P! q factor that drives the angular dependences in 
the former case there is a continuum of |p|, |q| values in 
the present case; for computational purposes we need to 
divide p space into a number finite size bins, just as we 
did the angular space in the simulations reported in [l4| . 
For our primary stability analysis we use only one bin, 
taking all neutrinos to be concentrated, say, at the peak 
of a thermal spectrum |p| = p Q . We have examined sim- 
ulations of the same problem with several bins and find 
no qualitative difference in the results for the instability 
conditions or results. 



1 Note that the spinor factor in the exchange amplitude is of 
the form [wi(p)it2(q]['U3(q)M4(p)]. We use a Fierz transform in 
reducing to the form shown in l|10[ l. We give the rather complex 
result, for the cases of interest, in sec. 6. 
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In the one-bin simulation the Hamiltonian is 

H cS = gnvNvXi £ ^S a Sp + gn v N u X 2 £ r^ p S a S p , 

(13) 

where 

Ai,2 = J dz P)q Di i2 (po,2>o,2;p,q) -1 (l - a; p , q ) , (14) 

and 

[S a ,Sp] = N- 1 ^2hl tP S y . (15) 

7 

4. N NEUTRINO CALCULATIONS 

In the mean-field (MF) reduction of equations ((Sj> for 
the Zq exchange problem, or the analogous equations for 
the scalar exchange case, derived from (fT3|) and (fT5|) . the 
particle number does not appear. We shall discuss so- 
lutions in the following sections. But insight is gained 
by solving the operator equations for systems with a fi- 
nite number of neutrinos. This can be done numerically 
for numbers N„ of a few hundred. One might have hoped 
that the rate corrections to the MF approximation would 
be of order iV^ 1 compared to the MF value. The situ- 
ation turns out to be somewhat more complex, but the 
calculations reported below do support the MFA in the 
domains in which we need it. Additionally, the finite N v 
calculations show fascinating connections between unsta- 
ble behaviors of the complete system with v oscillation 
terms turned off, and unstable behaviors in the MF ap- 
proximation when the oscillation terms are reinstated. 

These calculations are based on the observation that 
the collective operators, N„S±, N U S^ , which obey an- 
gular momentum rules, are sums of individual neutrino 
"spin" operators for the individual modes a. There are 
2N V basis states, but when the number of bins is small, 
the size of the problem can be much smaller. To be- 
gin with, the total "angular momentum" of each bin, 
S( Q ) . S^"), is conserved. Then if we have an initial state 
for the bin in which every "spin" is up, this "angular mo- 
mentum" has maximum value, characterized by quantum 
number j a . The basis of states for this one bin is thus 
effectively reduced from 2^' ' to N^. 

In the very simplest case Friedland and Lunardini j2(| 
(see also [27j) found an analytic solution in the limit of 
a large number of u's for a two-flavor model with no v 
oscillation terms. They take an initial state with two 
unidirectional beams of different flavors, at relative an- 
gle 9. The point is to calculate the time scale for the 
exchange of a macroscopic fraction of the flavors from 
one beam to the other. In the notation of the last sec- 
tion, we choose variables and S±, S3, T±, T3, where the 



S's and T's are the collective respective flavor matrices 
for the two beams. The operative neutral-current inter- 
action is then the form 

H = 2 3 / 2 GiV(S + T_ + S-T+ + ^S 3 T 3 ) , (16) 

where we have defined G = 2~ 1 / 2 Gftc„. In the standard 
model case described by © the parameter a is unity; this 
is the case discussed in ref. (26|. But the a < 1 case, un- 
physical here, is a simple prototype both for Z exchange 
models cases with complex angular distributions, and for 
all the results for scalar exchange that we shall describe 
later. 

We take the initial state to be an eigenstate of S 3 ,T 3 
with eigenvalues S3 = 1 and T 3 = —1; that is to say, the 
respective beams are of pure flavor and of different fla- 
vors. As displayed in ([7]) the equations of motion for these 
operators, when expresed in terms of neutrino densities, 
are free of any explicit occurrence of N v (or, equivalently, 
the volume). Nonetheless, in the products S^S^ on 
the right hand side of (J7J the operators act in the com- 
plete N„ neutrino space, and the solutions will have N v 
dependence. Solving numerically for the case a = 1, the 
mixing time is found to be of order, 

t ev ~ G~ 1 N^ 2 , (17) 

essentially the result of ref. 126] achieved by less elegant 
means. 

There can be no physics here, since this mixing time 
for a system of a very large number of neutrinos would be 
longer than that originating in the ordinary cross-section. 
Thus the truncation to the forward Hamiltonian, as one 
might have expected, threw away almost all of the actual 
flavor exchanges predicted by the full theory. 

In contrast, if —1 < a < 1, the results are very differ- 
ent. The analytic technique of Ref. 26] not being ap- 
plicable here, direct numerical calculation for these cases 
gives the mixing time, which we now explicitly define as 
the time at which £3 passes through the value zero, as 

h -CaG- 1 log N„ (18) 

The coefficient C a is of order unity for the fastest case 
a = 0. Eq. (Ti"5]) holds with somewhat smaller coefficient 
for \a\ < .9 in our numerical calculations. Again we are 
not at the moment suggesting a physical application of 
(p~8|) . though logarithms never really get large. 

In Fig.l we show results for values of ^=100,400,1600 
for the two cases a = .5 and a = 1. The logarithmic de- 
pendence on N v in (TT5)) is reflected in the equally spaced 
intercepts of the solid curves, and the \fN v dependence 
in (fi~7)) in the behavior of the dashed curves. 

We will not in this section try to analyze how these 
finite N results could reflect real behavior of neutrinos 
that are at large in space. In examples that we develop 
in the next section we find that, as one might even ex- 
pect, when the basic model is augmented with a tiny seed 
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FIG. 1: Flavor exchange between two beams as a function 
of time. The curves are labeled by the number of neutrinos 
used in the simulation. The dashed curves are for the case 
a — 1 and the solid curves are for the case a — .5. The time 
is measured in units G _1 . 

that can mix the states, the mixing time is given essen- 
tially by (HH|) but with log.ZV„ replaced by | log£| where £ 
measures the magnitude of the seed. For the case of Z 
exchange and flavor mixing, the seeds are the ordinary 
oscillation parameters. For the scalar exchange model 
and helicity reversal, we return later to the question of 
possible seeds. If nature provided none (except for some- 
thing so super-small that even the logarithm was huge), 
then possibly the log-/V„ model as is would be interest- 
ing. We would characterize the situation as "seeded by 
quantum fluctuations" . 



and the initial conditions 

S 3 (0) = 1 , T 3 (0) = -1 , T+(0) = T_(0) = , (21) 

the equations (|2"0)l determine the evolution of the system. 
In the absence of oscillations, e = 0, none of the variables 
changes in time for any value of the parameter a. 

In fig. 2 we plot the results of solutions of the set 
(|20| for the two cases, for three different values of the 
vacuum oscillation parameter e, with the range of the 
latter chosen to emphasize the similarities between fig. 2 
and fig.l. Extending the calculation to the entire region 




FIG. 2: As in fig.l, the flavor exchange against time, now 
from the mean field equations, and labeled by the strengths 
of the oscillation parameter e/G, plotted for a — .5 (solid), 
and a = 1 (dashed) 



5. MEAN FIELD CALCULATIONS 

Turning to the mean field approximation, we add an 
ordinary neutrino oscillation term H osc to the Hamilto- 
nian (fT5|) . 

ffosc = eN u (S+ + S-+ T+/2 + T_/2) (19) 

where we used different oscillation rates for the two dif- 
ferent bins, in recognition of the fact that in applications 
the oscillation rates depend on energy, and the bins will 
typically have different energies. We use the commuta- 
tion rules (JS|) among the S operators and their analogues 
among the T operators to obtain the equations of motion, 

i-TL S + = G(S 3 T + - aT 3 S + ) + eS 3 , 
at 

i^-T+ = G(T 3 S+ - aS 3 T + ) + eT 3 , 
at 

ij t S 3 = -2G(S-T+ - S+T-) + 2e{S+ - S-), 

i^-T 3 = -2G(TS+ - T+S-) + 2e(T+ - T_) . (20) 
at 

With the mean field identities, S- = S+, T_ = T£, 



10 6 < e < .04 we find a good fits to the values of the 
S3 = intercept t| v to be given for the a = .5 data by 

tT = G- 1 log[G/e] (22) 

and for the a — 1 data by, 

t™ = 1.85(Ge)- 5 , (23) 

In spite of the intriguing similarities between fig.l and 
fig. 2, we still at this point have learned nothing about 
whether or not the mean-field equations give the com- 
plete solution in the limit N v —¥ 00. Ideally, now, we 
would go back and add the vacuum oscillation terms to 
the calculations of the last section and try to replicate the 
curves of fig. 2. The technical problem is the following: 
for the case with e = and with our simple initial con- 
ditions there were only N v simultaneous coupled ODE's 
to be solved. When oscillations are included, there are 
equations. Furthermore, as we see from comparing 
the a = .5 plots in figs. 1, and 2, a value iV = 100 is 
not sufficient to push the transition region with e = 
out far enough to see the specific effect of e. Therefore 
our computational limit N = 100 for the case with the 
oscillation term included is insufficient for a test of the 
mean field approach for the values of e used in fig. 2. 
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However, we have been able to squeeze in a simulation 
for a value e = .2G, in the N v = 100 case, at times so 
small that the "finite N v " contamination hardly enters 
(and so short that ordinary oscillation for the G = case 
is negligible) and we find excellent agreement with the 
mean field result. We have also made a similar check of 
the MF approximation in a very limited domain in the 
a = 1 case. 

The key to the vast qualitative difference in the behav- 
iors for the a — 1 case and the a < 1 case for either of the 
above approaches is found by doing a standard instability 
analysis on the equations ([20)) . Setting e = and lineariz- 
ing in the perturbations S± , T± we see from the last two 
of the set that, to linear order, S3 and T3 remain at their 
initial values. Then the solutions for S±,T± have sim- 
ple exponential behavior e ±A * where A = — a 2 )GN, 
and the real exponential for the case a < 1 signifies an 
exponential growth in response to a perturbation, such 
as the e perturbation that gave rise to the rapid mixing 
in (f2"2"j). The slower mixing in the case a = 1 of (f2l?|) . is 
still much faster than the oscillation time e" 1 in the case 
that G/e = Gpii^/e >> 1, which holds in the regions of 
most interest in a supernova. In this case we can roughly 
say that two zero eigenvalues of the response matrix of 
the linearized system are effectively driven complex when 
e is introduced and the system evolves a little, leading to 
a speed-up but not as dramatic a speed-up. We char- 
acterize this behavior as "moderately fast" , and that of 
([22]) as "very fast" . 

In the application to neutrino physics in a supernova, 
there have been a number of numerical simulations 
llj that produce flavor exchange at a rate much faster 
than vacuum oscillation rates and occurring somewhat 
outside the neutrino-surface. These simulations must di- 
vide the neutrino flow into many bins, since the angular 
dependences coming from cos(0 p , q ) in ([3]) are critical. 
The authors do not analyze their results in terms of the 
stability matrix of the linearized problem, but it is our 
conclusion that the closer-in of these phenomena are of 
our category, "moderately fast", in the sense that they 
stem from a pair of zero eigenvalues. They also appear 
to fit the (G^n^e) -5 behavior of our simple case. 

However, there cases in which the "very fast" evolution 
appears not only possible but likely [l3j|. This can 

occur when we have multiple beams simulating an initial 
flavor dependent angular distribution more disorganized 
than those assumed in refs. [l[-[ll| (which have perfect 
cylindrical symmetry). This appears to be the case in 
the region just below the surface of last scattering for a 
typical v e . In studies reported in [14[ we started with a 



basic distribution in this region in which the and v T 
neutrinos (and anti-neutrinos) have angular distributions 
biased more outward that than those of v e , v e . With this 
distribution, we found "moderately fast" evolution, in a 
number of cases, with the most dramatic results for the 
"inverted hierarchy" assumption regarding the neutrino 



mass matrix. 

Then we added small random transverse distributions, 
breaking azimuthal symmetry. About 50% of the time 
the behavior shifted to that of very fast evolution, with 
mixing lengths much less than 1 cm. These results were 
robust with respect to introducing a neutrino energy 
spectrum, through the energy dependence of the oscil- 
lation terms, and the results were essentially the same 
for the normal and inverted hierarchies. In view of the 
turbulent fluid motion in this region, this suggested com- 
plete equalization of the energy spectra of the various 
neutrinos emerging from the neutrino-surface, with pos- 
sible effects on explosion dynamics and the R process, 
and with certain effects on the neutrino pulse received 
on earth. In contrast to standard MSW transformations, 
which are driven completely by the change in background 
electron-density changes, these instabilities only require 
time to pass in order to establish the prerequisites for 
later sudden changes, with all model parameters remain- 
ing time independent. 



6. SCALAR EXCHANGE MODELS 

With the results of the previous sections in mind we 
now investigate the behavior of the model described in 
section 3. Beginning from the one-bin equations (|13[) . 
(fT4|) and (p~5|) . we note first that in the mean field inter- 
pretation we can, for the purposes of deriving the equa- 
tions for the collective variables , cut down notation 
by eliminating all sums over modes p in their definitions, 
and at the same time setting N v — 1 in (IT31 and (|15l) . 
We give a symbol for each of the 16 bilinear forms in the 
operators a, 6, a, b, introduced in sec. 3, and then list 
all 120 independent relations of the commutator algebra. 
In this notation we have the 16 operators that we divide 
into three groups, 

Group I 

Qj, = b^b — a' a , Q3 = Pb — a' a , 
No = b^b + a f a + tfb + a^a , 
Z = ata + tfb-tfa-Pb. (24) 

Group II: 

Q+ = tfa , Q_ = a)b , Q+ = tfa , Q_ = a)b (25) 
Group III: 

X+ = tfa , X_ = a f 6 , Y+ = tfa , y_ = a)b . (26) 
Group IV 

V+ = a^a , V- = a) a , W+ = b r b , W- = b^b . (27) 



The Hamiltonian (|13|) is given in terms of these oper- 
ators by 

H cS = gn v Xi U(X+Y+ + X_F_) 
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+ (2 + \){V+W+ + V-W-) + A(A+A_ + Y+Y- 

Q+Q+ + Q-Q- - Q3Q3/2 + Q+Q+ + Q-Q- 
ji + \(Q+Q- + Q-Q+ - O3Q3/2) + iv 2 - vz 2 } , 



hQ 3 /2 + 

where A — Ai/A 2 , r\ — 1/4 
Ai = (2p )" 2 



A/8 and 
2 

log 



A 2 = (2po)" 



rn 

2 



4p 2 + TO 2 . 



- 1 



m 
ip 2 



lof 



r 4p 2 



(29) 



The equations of motion are given as usual by id/dt S — 
[S, H] where we now use the commutation rules that 
come directly from the definitions (|24p - (j2"7| . 2 

We shall always suppose that we start at t = from a 
state in which all of the operators that mix one helicity 
with another, or that mix neutrino with anti-neutrino, 
are zero. The reason for dividing the operators into 
the four groups is that there are two cases in which the 
problem simplifies a little: a) when the mean fields of 
all of the operators of group II and group IV , namely 
Q±,Q±,V±,W± remain zero; or b) when the operators 
of group III and group IV X± ,Y±,V±, W± , remain zero. 
In either case the linearized mean-field equations have a 
growing mode, when we take initial condition of all-active 
helicities. Both cases require a seed to make anything 
happen. In case b, we believe that a neutrino magnetic 
moment in a magnetic field can provide the seed, but 
discussing it in detail would require backtracking to put 
angle dependence in the equations and require binning in 
angle in the simulations. We treat case (a) in detail in 
the next section. 



7. THE CASE; Q± = Q± = V± = W± = 

As in the eariier mean-field examples, nothing will 
happen until we add an analogue of the neutrino os- 
cillation terms used above. In the present case we 
use the minimal lepton-number breaking bilinear, H 1 = 
j dxip(x)Cil>(x) ) where m v is the neutrino mass and 
C the charge conjugation matrix. This translates, in our 
notation, into a mixing term in the forward Hamiltonian, 



H' = e(A 4 



X. 



Y+ + F_ 



(30) 



This will catalyze reactions in which there is sudden to- 
tal exchange of helicity between the neutrinos and anti- 
neutrinos. Now defining G = gn v \\, the effective Hamil- 
tonian is, 



2 Before making the mean-field approximation, we could worry 
about the order of factors on the right hand side of the equations. 
But the commutator terms that distinguish one operator order 
from another would give corrections of relative order N^ 1 in the 
end, and thus there is no ambiguity in the mean-field equations. 



H = G 2X+Y+ + 2X-Y- + A(A + A_ + V+V_) 



(2#3Q 3 /2 - Q3Q3/2 - AQ3O3/2) - Z 2 



This, together with the commutation rules, 

[A + ,A_] = (Q 3 + Q 3 -Z)/2, 
[Y +1 Y_] = (Q 3 + Q 3 + Z)/2 , 
[X+,Q 3 ] = [X + ,Q 3 ] = -X+ , 

[y_,Q 3 ] = [F_,Q 3 ] = y_ , 

[X + , Z] = 2X+ ; \Y-,Z] = 2Y_ , 



H' 



(31) 



(32) 



determines the equations of evolution. We easily see that 
Q3 ~ Q3 is conserved, and we will always take initial con- 
ditions such that Q%(t) = Qz{t). With that substitution 
(EU and (EU) give, 



l dt X + = ° 



2Q 3 [V_ + (1 + A)A+] 



-Z[Y_ 



A, 



(4 ?7 +-)A + ] 



e(Q 3 - Z/2) 



i—Y^ = G 
dt 



2Q 3 [-A+-(1 + A)y_] 
-Z[X+ + (4ri + ~)Yl] - e(Q 3 + Z/2) , 



i—Q 3 = 2G[X+Y+ - A_Y_] 



e[X + - X_ + Y+ - Y_] , 
2A_ +2Y+ - 2Y~_) 



i— Z = e(-2X 
dt v 

As in the earlier examples, we can investigate the pos- 
sibility of very fast mixing by looking at the linear per- 
turbations away from the trivial solution with e = 0, 
X+(t) = Y-(t) = , and Q 3 (t) —const, by looking at 
the two equations, 



(33) 



i-X+ = 2GQ 3 [Y + + (l + \)X + ], 



ijY_ = -2GQ 3 [X + + (1 + A)Y1] 



(34) 



The solutions have time dependence exp[±wo^] where 
luq = 2GQ 3y / l - (1 + A) 2 . If -2 < A < we thus find 
real exponential behavior that can drive rapid mixing. 
Taking a value of A = —1 appropriate for most of our 
mass region and solving the full equations (|33|) we plot 
the results in fig. 3 for a set of values of the "oscillation" 
parameter e, ranging over eight decades. The intercepts 
with Q 3 = arc well fit by 



w o 1 l°g[ £ / w o] 



(35) 
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FIG. 3: Transformation of left handed v's and their an- 
tiparticles into right-handed counterparts, as catalyzed by a 
tiny lepton violating oscillation parameter e. The results for 
e/G — 10~ 3 are given by the heavy curve. Moving to the 
right, we show the plots for values, 10~ 5 , 10~ 7 and 10~ 9 . 
Time is measured in units of G _1 . 

Although we have lepton number violation, measured by 
Z, catalyzed by the term (|30[) . at no time in the evolution 
does Z become non-negligible, even for our largest value 
of e. 

Finally, we address the fact that the neutrino absolute 
momenta are to be distributed in a thermal spectrum, 
so that we must introduce multiple beams of neutrinos 
of different values |p|, exactly as we introduced multiple 
angles in the Z$ exchange model 3 . The effective Hamil- 
tonian is now , 



time evolution of J^. Q3 are the same as in the one 
bin model. But the scatter of values obtained in our 
procedure can be expected to have some effect. We have 
experimented in fitting the parameters X\ J for the case 
of an initial thermal distribution at temperature T with 
up to 9 energy bins, each containing the same number 
of neutrinos, then solving the equations derived from 
(|36| and comparing with the results of the single energy 
approximation, where we set the parameter po = 2T. 
There is never more than a factor of two difference in 
the turn-over time, and discrepancies of this amount 
are found only in the region I0~ 5 < m/T < 10 -4 . 
We have not attempted the calculation in the region 
m/T < 10- 5 . 

Fig. 4 shows an example of a multi-bin simulation in 
which the elements of the matrix vary over a factor 
of 6, as compared to a single bin simulation in which Ai 
is taken as an average. 



e 3 




- 1.0 p 



N B 
a,/3 



X 



(/3) 



,y(«)y(/3) _ 



W/2]) 



+e^[(x| Q) + X { _ a) ) ± (Yl a) + F M ) 



(36) 



Here the matrices Ai,A2 are derived from 
Jd[cos6'p,q]£)i i 2(p,q) where the functions D are 
from (fTTj) and $F2\i , evaluated on a mesh of points 
|p|j, |q|j. In the numerical simulation we have taken 
a subdivision into Nb bins such that with an initial 
thermal spectrum every bin is equally occupied. The 



=,(<) 



is then to be taken as 



initial value of each Q 3 , ^ 3 
1/Ng. If all the matrix elements of Aj 3 are equal to 
each other and likewise for A^, then the results for the 



3 Our angular averaging was legitimate, however; note that in the 
Zq exchange model we could have taken isotropic initial distribu- 
tions and averaged away the cos(0 Pi q) factor in the interaction. 
But there the interesting physics is for non-isotropic distributions 



FIG. 4: The quantities in the same units as in fig. 4. The 
dashed curve is for the single bin theory as plotted in fig. 4, 
here shown for a mixing parameter e/G = 10~ 4 . The solid 
curve is for the multi-bin model with coupling constants vary- 
ing over a factor of 6, rather than for a single G. 

In all of the situations discussed above, our universe 
contained only neutrinos and antineutrinos. We believe, 
however, that the results are applicable under some con- 
ditions to real situations in which there are other parti- 
cles that interact with neutrinos, such as the electrons 
and positrons in the early universe. A condition that we 
must demand is that the rate of neutrino reactions from 
ordinary scattering of neutrinos in the medium be small 
compared to the rate of transformation by the coherent 
mechanism. The limits on the scalar coupling that we 
shall quote demands that rate for the coherent process 
be at least five times that based on [(cross-section) x den- 
sity]. However, in the case of the scalar coupling model, 
there is a further issue that we must discuss before con- 
cluding that we do indeed have a model that predicts 
rapid creation of sterile helicity states. When we choose 
the scalar particle (0) mass to be substantially less than 
the temperature, which we took as 1 MeV in the appli- 
cation, we must worry about the real process v + v — > cf>; 
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for equal mixtures of active and sterile particles this re- 
action rate is exactly of the same order of magnitude as 
that of our rapid coherent process 

In the initial state in our calculation, occupied only 
by active neutrinos and anti-neutrinos, this process is 
forbidden, but as our state with mixed in sterile parti- 
cles develops, there would be reactions v + v s — > <j) and 
v+ v s — > (J). One would guess that this would not change 
the qualitative results, but it would have been good to 
have an analytic treatment of the effect. The issue, more- 
over, appears to be essentially the same as that posed 
in ordinary neutrino oscillation theory when incoherent 
scattering is introduced. This has been discussed many 
times, going back at least to ref. [3(J, and we can state 
a typical outcome as follows: in the case of a large an- 
gle neutrino oscillation where the neutrino scatters as 
well, with flavor conservation but flavor dependence in 
the scattering rate, the scattering makes little difference 
over a few oscillation periods if the two time scales are 
comparable. If the scattering rate is much faster than the 
oscillation rate than the oscillation rate in the absence of 
scattering, then the flavor mixing is slowed, even frozen 
(in the limit of large scattering). But in our analogous 
system, the incoherent process of <j> production is never 
strong enough to produce a major effect. 



8. DISCUSSION 

We have shown that the minimal Yukawa coupling 
of a Dirac neutrino to a scalar particle can, through 
a collective effect, lead to rapid production of other- 
wise nearly sterile right-handed neutrinos, in an ini- 
tial thermalized state of left-handed neutrinos and anti- 
neutrinos. The inverse time scale is essentially of order 
{g 2 T/ log[T~ 1 e], where e measures strength of lepton- 
nonconservation from some other source. For any rea- 
sonable value of e other than zero this is much faster 
than the rate of production g A T in an ordinary perturba- 
tive calculation of vl + Dl — ► vr + vr from cross-sections 
and density. 

If we extend the model to include all three neutrino 
species, universally coupled to the scalar, then we can 
use the estimate (|35|) to put limits on the parameters g, e 
to avoid complete population of the sterile (or "wrong- 
helicity") states, contributing an additional three neutri- 
nos to the mix in the period shortly before decoupling. 4 

The estimates do not depend critically on the mass of 



the scalar particle as long as it is in the range < m < T. 
In setting limits on g we demand that the characteristic 
time for the helicity overturning effect be less than 1/5 of 
the collision time for an ordinary reaction, like e + v — > 
e + v in the medium at a temperature of 1 MeV. This 
gets us to a limit of about g < 1 CP 10 , for values of the 
lepton number changing parameter e > 10~ eV. This 
limit is less than those that come directly from particle 
experiment, and from supernova analysis, the latter at 
the g < 1CT 4 level. 

Alternatively, if, as seems to be indicated by recent 
WMAP analyses [3], more that three light neutrino 
species are needed, then it would be easy to adjust pa- 
rameters to get any number between three and six in 
these models. 

If couplings of the scalar are taken to change neutrino 
flavor, then the stability analysis gets very complex but is 
still straightforward. However the purpose of the present 
paper was to point out some range of possibilities rather 
than to try to do detailed phenomenology, with the sec- 
ondary purpose to understand better the workings of the 
much discussed instabilities in the Z exchange case. In 
the latter case we have found that the mean-field assump- 
tion is supported by simulations involving a few hundred 
neutrinos. 

In our detailed model, on which the above estimates 
are based, the seed was taken as a tiny term that vio- 
lates lepton number conservation. It appears that a neu- 
trino magnetic moment in a local magnetic field could 
also drive the effect. Given various more fundamental 
parameter assumptions, of a kind practiced in the game 
of "model building" , one might be able to argue that one 
or the other should dominate. A third possibility, briefly 
discussed at the end of sec. 4 is that quantum fluctuations 
could serve as a seed. 

Other physical systems can show collective behaviors 
analogous to the ones that we have discussed for neutri- 
nos. Indeed an exactly analogous case has already been 
pointed out. Kotkin and Serbo showed 3l|, in a purely 
classical calculation, how a beam of polarized photons 
can exchange polarizations with another beam of polar- 
ized photons, through the (Heisenberg-Euler) quartic ef- 
fective [32| interaction in QED. This conclusion does not 
require that the beams be monochromatic or coherent, 
though in practice the intensities needed for an experi- 
ment would surely require laser sources. 5 As in the case 
of our very fast models, the calculated rate for the process 



4 Of course, in the model we used in this paper the creation of the 
sterile states depletes the occupancy of the active states. But in 
the actual environment, before v decoupling, the ordinary weak 
interaction rates, although slower than our transformation rate, 
will re-populate the active states; in the end we would expect 
rough thermal occupancy for all states. 



The required parameters do not appear to be very far beyond 
those of the SLAC experiment [32| that detected the reaction 
7 + 7 — > e+ + e_ , rather misleadingly reported as "photon- 
photon" scattering . Although few would doubt the correctness 
or applicability of the effective 4-7 coupling for photons in the 
ultraviolet, its experimental confirmation would be of some in- 
terest. 
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is many, many orders of magnitudes faster than the rate 
for polarization exchange as calculated from the text- 
book photon-photon scattering cross-section, by virtue 
of being lower order in the effective coupling constant. 
A derivation and generalization of these results follow- 
ing the approach of the current paper can be found in 
ref [33j]. This work was supported in part by NSF grant 
PHY-0455918. 
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